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Abstract: Moran’s / statistic, a popular measure of spatial autocorrela¬ 
tion, is revisited. The exact range of Moran’s I is given as a function of 
spatial weights matrix. We demonstrate that some spatial weights matri¬ 
ces lead the absolute value of upper (lower) bound larger than 1 and that 
others lead the lower bound larger than —0.5. Thus Moran’s I is unlike 
Pearson’s correlation coefficient. It is also pointed out that some spatial 
weights matrices do not allow Moran’s / to take positive values regardless 
of observations. An alternative measure with exact range [—1,1] is proposed 
through a monotone transformation of Moran’s I. 
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1. Introduction 

Spatial autocorrelation statistics measure and analyze the degree of dependency 
among observations in a geographic space. In this paper, we are interested in 
Moran’s (1950) I statistic, which seems the most popular measure of spatial 
autocorrelation. To the best of my knowledge, however, the property has not 
been thoroughly investigated. To be worse, there are some misunderstandings 
with respect to similarity to Pearson’s correlation coefficient. In some books as 
well as Wikipedia, it is claimed that the range of Moran’s / is exactly [—1,1] 
like Pearson’s correlation coefficient. However it is not true. 

In Section 2 of this paper, we give the exact range of Moran’s / as a function 
of spatial weights matrix. In our numerical experiment, some spatial weights 
matrices lead the absolute value of upper (lower) bound larger than 1. Others 
lead the lower bound larger than —0.5. Thus Moran’s I is quite different from 
Pearson’s correlation coefficient with exact range [—1,1]- It is also pointed out 
that some spatial weights matrices do not allow Moran’s I to take positive values 
regardless of observations. In Section 3, we propose an alternative measure with 
the exact range [—1,1], which is a monotone transformation of Moran’s I. 

Actually, this is not the first paper which points that the range of Moran’s 
/ is not [—1,1]. To the best of my knowledge, Cliff and Ord (1981), Goodchild 
(1988), de Jong, Sprenger and van Veen (1984) and Waller and Gotway (2004) 
did so. From mathematical viewpoint, Theorem 2.1 in Section 2 is essentially 
the same as the results by de Jong, Sprenger and van Veen (1984). However any 
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alternative with exact range [—1,1], which I hope, helps better understanding 
of spatial autocorrelation, has not yet been proposed. Furthermore it is notable 
that since the alternative is proposed by a monotone transformation of Moran’s 
I, the permutation test of spatial independence based on Moran’s I, which has 
been implemented in earlier studies, can be saved after the standard measure is 
replaced by the alternative. 


2. Moran’s I 


In this section, we give the exact range of Moran’s I. Suppose that, in n spatial 
units si,..., s n , we observe y\ = y(s i), ..., y n = y(s n ). We also suppose that 
spatial weights Wij, the weights between each pair of spatial units s* and Sj, 
which satisfy 

> 0 for any i ^ j, wu = 0 for any i, (2.1) 

are given by some preset rules. Then Moran’s I statistic is given by 

T _ n EHi E”=i w ij(vi - V)(yj - V) 

{EH=i^}EHi {yi-y ) 2 ' 

Let y = (yi,... , y n ) T and y = y - yl n = (/„ - l„l^/n)y where T denotes the 
transpose and l n denotes the n-dimensional vector of ones. With the spatial 
weights matrix W = {Wjj ), Moran’s I statistic is rewritten as 

_ n rWy 

E,,/ : "u v T v 

When W is not symmetric, the symmetric alternative (W + W T )/2 can be 
identified to W since 


y T Wy = y T {(W + W T )/2}y } for any y, 

E” ■ , = lT n W1 » = lT n{( W + ^ T )/2}ln = , {w ZJ + Wji }/2. 

Thus, in this paper, we assume that W is symmetric. 

Let V be the (n — 1)-dimensional subspace of R n orthogonal to l n . Then 
I n — l„l^/n is the matrix for projection onto V. The spectral decomposition of 
I n - l n ln/n is given by 

I n - 1„1 l/n = hihj = HH T ( 2 . 2 ) 

where hi € R ra , hjhi = 1 and hjl n = 0 for any 1 < i < n — 1. Further, in (2.2), 
H = (h \,..., h n _ i) is the n x (n — 1) matrix where {hi ,..., i} span the 
subspace V as the orthonormal bases. 

Let v = H T y G R n 1 . Then Moran’s I statistic is rewritten as 


y T Wy 

nwy T y 


y'HH'W HH'y 


v T Wv 


nwy T HH T y 


v T v 


(2.3) 
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where w = Ylij =1 w ij/ n2 and W is the (n — 1 ) x (n — 1 ) matrix given by 


W = (nw)~ 1 H T WH. (2.4) 

For (2.3), the standard result of maximization and minimization of quadratic 
form is available. Let Am and A( n _i) be the smallest and largest eigenvalue, 
among n — 1 eigenvalues of W, respectively. Then, for any v £ K n_1 , 

I € [A(!), A(„_1)] . 

The lower and upper bound of I is attained by 


j = J A (i) v « i(i) 

|^ A (n— 1 ) V (X j( n — 1)5 

where j(-|) and jt n -i) are the eigenvectors corresponding to A^) and A( ra _i) 
respectively. By the following equivalence, 

vocj <S=> Hy oc j ^ H Hy x Hj <^> (J — l n lf/n)y oc Hj 
<t=> y = dl n + c^Hj for any C\ and nonzero C 2 , 


we have a following result. 

Theorem 2.1. The exact range of Moran’s I is 

I € A( n _!)] 


where A^) and A( n _i) be the smallest and largest eigenvalue of among n — 1 
eigenvalues of W given hy (2.4). The lower and upper bound is attained as 
follows; 


A(i) y = Cil„ + c 2 iTj(i) 
A(n— 1) V = C3I n T ip 


where any Ci,C 3 and any c 2 ,C 4 7 ^ 0 . 


Remark 2.1. In (2.2), the choice of H or hi ,..., is arbitrary. Clearly A(i) 
and A( ra _i) are also expressed by 


A(i) = 


. z T Wz 

mm —--, 

z T i„=o nwz T z 


A(n— 1) 


z T Wz 

max —-- 

z t i r =o nwz T z 


which do depend on W not on H. In practice, when A( X ) and A( ra _i) are calcu¬ 
lated, the choice of H = (hi ,..., h n -i) based on Helmert orthogonal matrix is 
available, that is, hi for i = 1 ,..., n — 1 given by 

h l = I 1 ?,-*. 0 l_i_i)/iji(i + l). 
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Table 1 

attainable lower & upper bound of Moran’s I 


n\q 


1 


2 


3 


lower upper 


lower upper 


lower upper 


10 -1.066 0.935 

20 -1.041 1.006 

30 -1.029 1.013 

40 -1.023 1.014 

50 -1.018 1.013 


-0.541 0.831 

-0.526 0.981 

-0.519 1.005 

-0.514 1.011 

-0.512 1.012 


-0.482 0.746 

-0.457 0.955 

-0.449 0.995 

-0.444 1.006 

-0.441 1.010 


Example 2.1. For simplicity, suppose si,..., s n are on a line with equal spacing 
and spatial weights for .s l: Sj is 



Table 1 gives lower and upper bound of Moran’s I for some n and q. For any n 
and q , the range of I is not [—1, 1 ]. It is demonstrated that some spatial weights 
matrices lead the absolute value of upper (lower) bound larger than 1 and 
that others lead the lower bound larger than —0.5, which is unexpected. Thus 
Moran’s I is unlike Pearson’s correlation coefficient with exact range [—1,1]. 

Remark 2.2. In all cases in Table 1, both A(!) < 0 and A(„_i) > 0 are satisfied. 
Researchers and practitioners would like measure of spatial correlation to take 
both positive and negative values depending upon positive and negative spatial 
correlation. It is clear that indefiniteness of IV is equivalent to A( X ) < 0 and 
A( n _i) > 0 simultaneously. Actually, as in below, some spatial weights matrices 
lead A( n _!) < 0, or negative definiteness of W. In such a case, Moran’s I cannot 
take positive values for any y , which is not desirable at all. 

In order to investigate definiteness of W, the following properties are useful 
on eigenvalues, trace or sum of all diagonal elements and their relationship; 

PI When AB and BA are square matrices, tr(AB) = tr (BA), 

P2 When both A and B are square matrices, tr {A + B) = tr A + tr B, 

P3 Trace is equal to sum of all eigenvalues. 

From PI, P2 and ^ wa = 0, we have 


trIP = (nw)~hrH T WH = ( nw)-hiWHH T 

= (nw)~ 1 trW (I — 1„1 n/n) = (nw)~ 1 trW — (n 2 w}) - 1 tiT)(TPl„ ( 2 . 5 ) 



From P3 together with trlP = —1, A^) < 0 are guaranteed whereas A( ra _!) > 0 
are not. 

In fact, if Wij = 1 for any i 7 ^ j or W = l n l^ — I, we have nw = n — 1 and 


W = (nw)~ 1 H T WH = (n - l)' 1 fl' T (l„l^ - I)H 
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= -(n - l^H^H = -(n - l)~ 1 I n -\ 

whose eigenvalues are all equal to — (n — l) -1 . This implies that W is negative 
definite and that Moran’s I is negative for any y. 

The negative definiteness with Wij = 1 for any i ^ j suggests possible negative 
definiteness for with small variability. Suppose, for 0 < a < 1, we observe 
n{n — 1 ) uniform random variables on (1 — a, 1 + a), which are assigned to Wij 
for i ^ j. We are interested in definiteness of 

W = {nw)~ 1 H T {(W + W T )/ 2} H. 

From numerical experiment, we have negative definiteness and indefiniteness of 
W for a £ (0, a*) and a £ (a*, 1 ), respectively, where a* is approximately equal 
to 0.3 for n = 25; 0.2 for n = 50; 0.14 for n = 75; 0.12 for n = 100. 


3. An alternative to Moran’s I 


In section 2, we pointed out some disadvantages of the original Moran’s I. In 
this section, we propose an alternative to Moran’s I, the range of which is 
exactly [—1,1] and which is a monotone transformation of the original Moran’s 
I. There are two steps of modification. The one is a linear transformation for 
acquirement of indefiniteness. The other is an adequate normalization for the 
linear transformation. 

As pointed out in (2.5), the sum of eigenvalues of W = (nw) -1 H T WH is 
— 1, which seems related to the bias of the original Moran’s /, that is, 

m = - ~~r < 0 

n — 1 

under spatial independence, as in Cliff and Ord (1981). In order to get a kind 
of unbiasedness, we start with a linear transformation of Moran’s J, 


(n- 1)1 + 1 


(n 


1 ) 


v T Wv 

V T V 


+ 1 = 


v T Qv 

V T V ’ 


where v = H T y £ R n 1 and Q = {n — 1)FF + I n -\. Since we have 
trQ = (n — l)trW' + tr/ n _i = —(n — 1) + (n — 1) = 0, 


the (n — 1) x (n — 1) matrix Q is indefinite, that is, (n — 1)1 + 1 can take both 
positive and negative values depending on y. 

We are now in position to give the adequate normalization of (n —1)1+1. The 
smallest and largest eigenvalue of Q is (n—1)A(!) + 1 and (n— l)A( n _!) + 1, which 
are guaranteed to be negative and positive respectively, by the indefiniteness of 
Q or trQ = 0. Hence, for v which satisfies v T Qv < 0 or equivalently (n— 1)1 + 
1 < 0 , we have 



l)A(i) +1 


v T v 
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and the equality is attained by 7 = For v which satisfies v T Qv > 0 or 
equivalently (n — 1)7 + 1 > 0, we have 

v T Qv , . , 

—~— < (n— 1)A(„_ 1 ) + 1 

v T v 

and the equality is attained by I = A 

Let Jm be an alternative to Moran’s I given by 

I _ (n — 1)7 + 1 _ f \(n — 1 )A(r) + 1 | (n— 1)7 + 1 <0 

C ’ + 1 (n—1)7 + 1 >0 

where An) and A^.j) are the smallest and largest eigenvalue of the (n — 1) x 
('n — 1) matrix W given by (2.4). Then we have a following result. 

Theorem 3.1. The exact range of Im is 


Im € [—1,1] . 


The lower and upper bound is attained as follows; 


-1 I = A 
1 7 = A( n _i) 


As far as we know, this is the first measure with exact range [—1, 1 ], which 
I hope, helps better understanding of spatial autocorrelation. Furthermore it is 
notable that since the alternative Im is a monotone transformation of Moran’s 
7, the permutation test of spatial independence based on Moran’s 7, which has 
been implemented in earlier studies, can be saved after the standard measure is 
replaced by Im- 
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